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Abstract: 

Tungiasis refers to the infestation caused by the permanent penetration of the female sand 
flea "Tunga penetrans" into the skin of the human or animal hosts and causes cutaneous 
lesion. In this paper, we developed a deterministic model for the dynamics of tungiasis in 
communities of human beings, animal reservoirs and flea infested environment in order to 
understand the way tungiasis disease spreads in a poor resource community. We 
established the conditions for local and global stability of the disease-free and endemic 
equilibrium points. The computational results showed that the disease-free equilibrium 
point is locally and globally asymptotically stable when the model basic reproduction 
number R, is less than unit, that is R, «1. Using the Lyapunov stability theory and LaSalle’s 


Invariant Principle we found that the endemic equilibrium point (EE) is globally 
asymptotically stable when R, »1. The numerical simulations have been presented to 


illustrate the way the model variables behave when there is no intervention and that the 
endemic equilibria exist and they are stable. On-host and off-host control measures which 
include focal spraying of insecticides to the premises, application of insecticidal dusts on 
animal bodies, environmental and personal hygiene, educational campaign, use of closed 
foot wear and application of plant based repellents should be implemented. 


1 Introduction 

Tungiasis is an ectoparasitic disease caused by the permanent penetration of the female 
sand flea "Tunga penetrans" into the epidermis of its host. Besides humans, various 
domestic and sylvatic animals are affected (Heukelbach and Feldmeier, 2002). Today, 
tungiasis is endemic in many countries in Latin America, the Caribbean and sub-Saharan 
Africa (Heukelbach etal. 2001). Transmission occurs when skin of the host comes into 
contact with soil or floor where adult female sand fleas have developed or when the host 
lives in close contact with infested animals (Pilger et a/, 2008). Models lead to a better 
understanding of the dynamics of the epidemic and the magnitude of the problem (Kaplan 
and Brandeau, 1994). Mathematical models can also help in figuring out decisions that are 
of significance importance on the outcomes and provides comprehensive examinations that 
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enter into decisions in a way that human reasoning and debate cannot (Bawa, 2013). 
Generally the disease dynamics requires a variety of mathematical tools, from model 
creation to solving differential equations to statistical analysis (Pathak et al, 2010). 


The basic concept in mathematical modeling is stability analysis of the equilibria of the 
epidemic models. The method based on the use of stable Metzler matrices has been 
proposed and proven to be useful to establish the global stability of the disease free 
equilibrium as described by Kamgang and Sallet (2008). The importance of the Metzler 
matrices is well recognized in the stability of dynamical systems and positive systems, 
more generally in biology, engineering and economics (Gupta, et al, 2014). On the other 
hand Lyapunov Direct Method by Korobeinikov (2002, 2004, and 2007) and McCluskey, 
(2006) combined with LaSalle’s Invariance Principle (LaSalle, 1976) has been a powerful 
tool for the analysis of stability of autonomous systems of differential equations through 
construction of suitable Lyapunov functions. Lyapunov functions are developed to establish 
the global asymptotic stability of the endemic equilibrium points. 


In this paper we consider the dynamics of tungiasis that involves the interaction between 
humans, animal reservoirs and flea contaminated environment. We determine the model 
equilibria and analyze their stability at local and global level. We establish the local stability 
of the disease-free equilibrium point using trace and determinant criterion, the global 
stability of the disease-free equilibrium point using stable Metzler matrix theory and we 
also establish the global stability of the endemic equilibrium point using Lyapunov direct 
Method combined with LaSalle’s Invariance Principle. Finally, we perform numerical 
simulations of the model system in a closed population to present the results. The rest of 
this paper is organized as follows. In the second section, we present the model with its 
basic properties. In the third section we carry out a qualitative analysis of the model 
whereby stability conditions for the disease-free equilibrium and the endemic equilibrium 
are derived. The fourth section presents different computer simulations of the system. In 
the last section, the biological significance of our analytical and numerical findings is 
discussed. 


2 Materials and methods 
2.1 Model formulation 
The total human population at any timet, denoted by М, is subdivided into three distinct 


epidemiological subpopulations namely susceptible humans S„, people who are mildly 
infested by jiggers /,, and the people who are severely infested by jiggers/,,,. 5, is 


generated through birth at a rate b, so that the human recruitment is b, N,. S, may 
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acquire infestation from the infested animal reservoir /,, and the flea infested 
environment F, following the force of infestation ¥,,and move to either /,, or /,, at the 
rates Pay 1,,/ М, and ag, Beute Е, /(К + F,,) respectively. Г, may acquire infestation from 
the environment as well and progresses into /,,, ata rate oz, P gyfr Fg /(k + Fe): We assume 
that individuals in classes S„ and /,, suffer only natural death at a rate „/„ and for class 
І,, the individuals suffer both the natural death at a rate „/„ and disease induced death at 
a rateo,,. Similarly the animal reservoir population denoted by М, is subdivided into 
three subpopulations namely, the susceptible animals reservoirs S,, animals that are 
mildly infested by jiggers J, and animals that are severely infested by jiggers/,,. 5, is 
generated through birth at a rate b, so that its recruitment rate is b,N,. S, may acquire 
infestation from infested animal reservoir /,, and the flea infested environment F, 
following the force of infestation ¥,, and move to either /,, or 1,, at the rates 
Pala, [Ng and ag Bere Е, /(k+F,) respectively. /,, may acquire infestation from the 
environment as well and progresses into /,, at a rate d, ғ, Ё, /(k + Еһ). We assume that 
individuals in classes 8, and 1, suffer only natural death at a rate wz, and the individuals 
class in 1,, suffers both the natural death at a rate „4, and disease induced death at a 
rateo,. The sub-model of environmental component consists of larvae and adult fleas 
compartments denoted by L, and Р, respectively. From within bodies of /,,, and 1,, 
eggs are shed by the gravid female flea into the environment at an average rate Ó, with 
contribution rates of 6,1,, and 6,1,, respectively so that they hatch and mature into adult 
fleas. The compartment L, suffers a natural mortality rate и, and matures to adult fleas 
class F, at a rate у,. At class F, the fleas are recruited from larvae as they mature and 
from the severely infested animal reservoirs /,,as they are shed at a rate e, with 
contribution rate €,/,, to the soil environment. The fleas are removed from the 
environment F, at the rate r, F; /(k + F.) and w, for burrowing into the host and natural 


death respectively. The forces of infestation for human population is given by 
Vu = Pag Lan Мн +U eu Beute Е, /(К + F,) and that of animal reservoir population is given 


by Y, =p, Та Бағы F, (k  F,). The variables and parameters of the model are 


summarized in Tables 1 and 2 respectively. 
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Table 1: The state variables of the model 






































Variable Description 

S, (t) Number of humans in a susceptible class at time, t 

T(t) Number of humans in a mildly infested class at time, t 

It) Number of humans in a severely infested class at time, 
t 

5 (t) Number of animals in a susceptible class at time, t 

I „ (t) Number of animals in a mildly infested class at time, t 

Lt) Number of animals in a severely infested class at time, t 

F, (t) The density of jigger fleas in the environment at time, t 

L, (t) The density of jigger larvae in the environment at time, 
t 

N „ (t) Total human population at time, t 

N (t) Total animal population at time, t 








Table 2: The parameters of the model 














Parameter Description 

k, K Half saturation constant of the jigger fleas and Environment 
carrying capacity of jigger larvae 

V ï Maturation rate from larvae to adult jigger fleas 

оор disease induced mortality rates for humans and animal reservoirs 
respectively 





Hg Hais Ups i, | Natural mortality rates for humans, animal reservoir, jigger fleas 
and jigger larvae respectively 





Vp The rate at which the jigger fleas leave the soil to attack the hosts 





Bru» Pea Effective contact rate between contaminated environment and 
susceptible humans, Effective contact rate between contaminated 
environment and susceptible animal reservoirs respectively 





Pau» Pa Effective contact rate between animals with fleas and susceptible 
humans, Effective contact rate between animals with fleas and 
susceptible animals respectively 




















bu ba Recruitment rates for humans and animal reservoirs respectively 

E, The rate of jigger fleas contribution by the severely infested 
animal reservoirs into the environment 

б, The rate of deposit of jigger eggs into the environment 

Á pr +: Ша The proportions of jigger fleas that leaves the environment to 
infest the susceptible human and animal reservoir hosts 
respectively 

М, The forces of infestation for humans and animal populations 








respectively 
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2.2 Model flow chart 

Using the above assumptions, definition of variables and parameters, the model flow 
diagram that depicts the dynamics of tungiasis transmission for the human population, 
animal reservoir population and the jigger fleas in the environment is shown in Figure 1. 


Animal гезеглог population 
b, N A 


Humans population 
b, Ny 


Jigger flea population 





(4t; + Fy Ms be 
\ 
\ 
\ 
KEY: \ 
> Flow 
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Figure 1: Tungiasis basic dynamical model 


Figure 1 shows the possible interactions between humans, animal reservoirs and jigger 
fleas in the environment. We have the susceptible human S4, infested humans at mildly 


and severe states /,, апа Г, susceptible animal S,, infested animals at mildly and severe 


states / ,, and/,,. We have larval class L, and the adult jigger flea class F}. 


2.3 Model differential equations 
From the compartmental model Figure 1 the following dynamical systems are derived to 
describe the dynamics of the transmission of Tungiasis: 
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Dynamicsin human population 
dS I 
=> =byuN y - [ou м; + Onn Brule ls, - My Sy 
dl I 
25 = Pan N 59 - Ad pu Beute K+ F, Lw - H yi Ï m 
dlm = QU pu P entr - E E Tu +U gu Beu Tr - È S, МЖ (ин +6 Mas 
dt EE. k+ F; 
Dynamicsin animal population 
dS I F 
25 =b,N, [od Ë), - H S a (1) 
dI Ts 
2s = 0, N, $4 — OG ga Pap k+ F lu — HaT a 
dl F F 
7 = Q rap gale k+ F ^ + Gea P gary ber" Е (ш, *O, Dm 
Dynamicsin Jigger flea population 
dL L 
P7 -д, aa s L,)- 6 + 4; Lp 
dF F 
та =J Ly t+ Ely, Mp Fg "r E+ F, 
where 
Ny(t)  Su() +10) - HF), NzO=S,041,,041,,0; 
F A F 
V gy = O gu Beute Lu Wan = Pan N V EN = Q ga D gale ERES 
D, 


Ил = Pay, Cn + Og = 1, 0<а;, <1 and 0<a,, <1 
A 


with initial conditions 


S,, (0) > 0,1, (0) > 0,9, (0) > 0, S , (0) > 0,7,,(0) > 0,2, (0) > 0, L; (0) > 0, F (0) > 0 


3 Properties of the model 

3.1 Invariant region 

To test whether the model is well posed epidemiologically and mathematically, we 
investigate the feasibility of the model solution whereby we present the model system (1) 
in compact form as in equation (2). We consider the approach used by (Mpeshe et al, 
2014). 
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m A(X)X +F (2) 
dt 
where X = (S4, I, Lg Sa Las lan Le Fe) > (3) 
A(X ) is the8 by 8 matrix : 
-D 0 0 0 0 0 0 0 
A, -D, 0 0 0 0 0 0 
A, A, -D 0 0 0 0 0 
0 0 0 -D, 0 0 0 0 
A(X) = (4) 
0 0 0 A, -D O0 0 0 
0 0 0 A A, -D 0 0 
0 0 6, 0 0 6 -D, 0 
0 0 0 0 0 в, Y -D 
where 


I F F 
D, = Pay `Ê +A pu Beute -_ E — Lg D, = Q py Вен! р; e *c,) 





Ny ІСЕЕ, k+ F; 
I F F F 
D, = Рау ыйыт? pep, Ies = APB rate per, Hes =Q ga Bea ГЕН 
д, r 1, І 
D, Èm +1) + ка) = Up oR = Pg N, ‚А, =Pa y 


F F F 
A, = а, Beute pp = Gy Due pU = Gg afr pir = (и, + c) 


and F isthe column vector such that: 

Е = (b, N, ,0,0,b,N,,0,0,0,0)" 20 (5) 
A(X) is a Metzler matrix for all X e RË considering the fact that F 2 0, model system (1) 
is positively invariant in RË, and F is Lipschitz continuous. Thus the feasible region Ф for 


the model system is the set 


Жиз о а аа s lal JOR : Sy S Nu Ig S Мк, 
D =J SN Í Sa SNL ЛА SN Sg HI r SN (6) 
SATI EN EKIPE LE 
Hr 
Therefore it can be verified that the solution remains in the feasible region Ф if it starts in 
this region. Hence, it is sufficient to study the dynamics of the model in D. 


8 STABILITY ANALYSIS OF THE DYNAMICS OF TUNGIASIS TRANSMISSION 


3.2 Positivity of the solution 
Let the initial data be ((S,, (2), /,„ (0.1, (0. S 4 (0.1 40.1 4,0.L-(, F„(t))>0) e & then the 


solution set ÍS 4 (t), Din (0),1,,(6),8,(0),1,,(0),1,,(0,1.:(0,Е0)) of the model system (1) is 


non-negative for all t > 0 


From the first equation of the model system (1) we have 


dS Ti 
"3 =bu М, 42, зу + Фен іні ls, - Hu S u 
H E 
dS I F 
En > Lou ааа p. a ys (7) 


Integrating equation (7) by separating the variables we have; 


' dS y Р Iu Fz 
Ta E A [ou N, + Onn Bene k4+F, + Ly, ds 


(01% F 
2 (orte AE tug E 
" ИДЕТ! 


S, (t) > S „ (Ое 
By using the similar procedure, it can be shown that the remaining variables 
Tuo 1н» Salan 1л Lg, Fgare also non-negative for all time т >0. Thus the solution set 
Ф = (5, (0), Lig 0, L4, O, S O, Lu), Lu (),L„(t),F„(t)) of the model system (1) is non- 


negative for allt > 0. 


4 Model analysis 
4.1 Existence of disease free equilibrium point 
The disease free equilibrium (DFE) point of the model system (1) is given by 


ЖЕЛ 





0,0, 


H HA 


Qe I soa alos) - 4 оодо, 


where by №, and № are constants with estimated values of 1500 humans and 1200 


animal reservoirs respectively. 


4.2 Local stability of the disease free equilibrium 

The local stability of the disease free equilibrium point is investigated using eigenvalues or 
trace-determinant criteria of the Jacobian matrix which is defined as a matrix of all first- 
order partial derivative of a vector-valued function (Simon and Blume, 1994). An 
equilibrium point is locally asymptotically stable if the Jacobian matrix evaluated at that 
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point has a negative trace and a positive determinant or has negative eigenvalues (Mpeshe 


et al, 2009). We prove the following Theorem 


Theorem 1: The disease free equilibrium ( ` ) whenever it exists is locally 
asymptotically stable if R, <1 and unstable otherwise. 


Proof 


At disease free equilibrium point D” we obtain the Jacobian matrix as given in (8) as 





- Uy 0 0 0 0. å -J 0 HA 
O -uy 0 0 0 J, 0 0 
0 O -J 0 0 0 0 Js 
д") 0 0 O -u, 0 -J, O -Jo 
0 0 0 O -u, J; 0 0 
0 0 0 0 D. m 0 Ji 
0 0 Ô, 0 0 à, -J, 0 
| 0 0 0 0 0 5, УУ, -Jo 
where 








(8) 


b b 
J, = (и, +01 )J, = Ј, = н „Ја =J; e А Je =(u,+0,)} J; = (у, +), 


H A 





J =J _ Ary Beute Nyby J =J Brake N aba J x le 
8 9 Кин ^" 10 11 ku, »* 12 F k 


The trace of ЛФ") is negative and its determinant is positive. Therefore all its eigenvalues 


have negative real parts. 


Tr (o )\< 0 which in this case implies 


-bu, +O, EM, tout ди, e 0 


Der Дә” ) >0 which in this case implies 


—4 Ш 
k 5, {Cen ен N uby ША (ш, +O, )+ «клк aba (uy +O7 
d gs Pg, N Ab €, | Vr 


| j о a y 





where R, = 


| is our basic reproduction number. 
pan + =) Mpk + Vp 


(9) 


(10) 
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We let; 
Ò = My (Un +O y Xr, + LU, us +0, Xu +r, X1=R,) 
О=6,уџ {Жкн Вън N ubu На (ш, TOS )+ Og Pg, N Рд (ин +O y ) 


For the determinant of the Jacobian matrix JÍ 2 to be positive © must be greater than О, 
that is @>Q and the basic reproduction number R must be less than one (R, « 1). It 
follows that the condition Detlj( °)>0 implies that TrU( ЕБ 0. Therefore, the disease 
free equilibrium (Фе) is locally asymptotically stable if and only if R, <land thus we һауе 


proved Theorem 1. 


4.3 Global stability of the disease free equilibrium 
The global stability of the disease free equilibrium is determined using Metzler matrix 
stability method by Castillo-Chavez et a/ (2002) whereby the system is put in the form: 








ты ӘК Ж 5:50 (11) 
| | 
dX. 

В.Х. 12 
dt 2 i ( ) 


where X, and X, are vectors for non-transmitting and transmitting classes 


respectively, X is a disease free equilibrium, and B,, B, and B, are matrices. For global 


DFE,n B, ? 


stability of the disease free equilibrium to hold the matrix B should have the eigenvalues 
with negative real parts and matrix B,need to be a Metzler matrix (i.e. the off-diagonal 


elements of B, are non-negative, denoted by B,(X)20,Viz j). 


From model system (1) we consider the following sub-systems 
X es S y. X, EU ызны ЫЕ; у. 


bu Ny b, N, (13) 
X DFEn (un No} 
Hg HA 


We consider the system (11) and differentiate it with respect to S,,,S, and L, at disease 
free equilibrium point Ф° we have the matrix B,,with which the eigenvalues are of 


negative real parts as indicated in equation (14) 
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“Ин 0 0 
B= 0 - 0 (14) 
0 0 -(у,%ш) 


Differentiating system (11) with respect to 1,,1,,,1,,1,,,Е; we obtain matrix B, іп (15) 








0 0 0 Pann d gu ВкнТЕЁРн № 
Hg Кин 
B,=|0 0 0 Paba Ga Be TE DAN А (15) 
Ha ku, 
0 6, 0 б, 0 





L 


Again from system (12) we differentiate it with respect to /,,,1,,.1,,,1,,, F, we have 


matrix B, as indicated in equation (16) 








=. 7 0 0 Pann 0 
Hg 
Он D gy DIN. 
0 РЕ +6 0 0 ЕНГЕН" F"H'H 
(ши " m 
b 
B,-| 0 0 -д, P 0 (16) 
A 
0 0 0 -(u "m ) d gj Bg Te DAN А 
a TOA 
ku, 
0 0 0 с, e 





We have investigated that the eigenvalues of matrix B, are real and negative and the 
matrix B, is a Metzler matrix therefore the model system (1) is globally asymptotically 


stable which is supported by the following theorem. 


Theorem 2: The disease-free equilibrium point (0°) is globally asymptotically stable if 
R, <1 and unstable if R, >1 


4.4 Existence of endemic equilibrium points 
Endemic equilibrium point Ф” is a steady state solution in which the disease persists in the 
population i.e. Zm = 0,1, £0,7,, £0,I,, #0and F, #0. We equate all expressions for the 


equations in the model system (1) to zero and let Ф” = (57, Г, Г. Sh Гы, Tan Le Е) and 
P =r, F; /(k +F) considering that the animal population is constant for b, > u, then 


from model system (1) at arbitrary equilibrium solution Ф we have. 
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P' =r, Fz /(k+F;) (17) 
E м.а) (18) 
А 
2 

Si = o bu (М,) = (19) 

Pau NA (b, — HA )+ СА ен ВенР + Un 
L = Peu N b, — Ha ba (N p Y (20) 

Nu (а, BP + Uy Voa NA (b, = H,)+ С. ен Ben P + c uy) 

N (b, — 
а, Brn Pb, (М, y te | к | 
D E +O,Ny EX + 2: Q1) 
4 PauN lb, H 
(ин +O H Wn (скн Ben P a 4) iol. : | 
“С, ен ВР +O 4Hy 
2 

5, = o b, (N.) - (22) 

РАУ, (b, — Ha )+ O JO paBriP t O Ag 
f o PsN b. (b, - Ha) Q3) 

Al = = 
lo. ou BP +O Ma Nex (b, — Ша) t+ Oye ЗР Og ) 
52 " 

jp E (24) 

QP «RP +T 
where 


D- б,КМ 4o («кн Ben Y (N.G, - Ua) Uy toy) +b„N2) 
U =5 KN N (ba — ua tg +O y (Pan N (ba -u,.)+ 20 My) 
=Ò KN yU gn Бен 

+b,Ny (сан NA (b, = 4)* uy) 
N (ba — Ma) My + On (олн (b, F u,)* 20, Hy ) 
+byNy (сан N (ba Е Ha)+ Hy) 

M = SKN „N , (b, “Ша ie (ин tOy Xo. N (ba - u,)) 

Q- (den Ben Y N„o, (6,b„N2 tug +O y )(6,N , (b, Е) И ш.) 

9, pau N (ba — Ha XN, G, - ua) uy t Og) шн Nu) 
R=N yX gu Вен w Е (us *Og Xv, * ш Холна (Ьл 7 шл)+ O AHg `) 
e 26,uy М, (b, т, Ha) 

Т-М,ш, (Us +O H Хода (ba - u, )+ Only NEN n (b, - Hi) +o +и,)) 

Substituting P` and Г,, from equations (17) and (18) respectively into equation (25) to 


U= Ô KN yQ gu Вен 








solve for F, such that 
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ж ж ж F. 
Oo Yale е Hele le T f (25) 
F; = lo 7 Ee +E, N (b -u,)-0,P | (26) 
O Arg 


Substituting F; from equation (26) into equation (17) we have 


* On Ye + тє N (ba = шл) слтР" 


P= (27) 





сш абу 2 +€,N (b, _-u,)-0,P' 
BP +0,y, LP" +AP -o0,P = суть +r,A- 0r, P 
(c,y,P^ -C)L, = A+o,P? - BF (28) 
where A= rue N (b, —ил) 
В - c uk (Б, = ua )* СА 
С= слу, 


Substituting L, from equation (24) into (28) we have 


*2 * 
бр dq +UP M) лаура 





QP” + ЕР «T 
loy P -c)pr? HUP +M)= (Ao, P? — BP' Jor" + RP +7) 
o QP” «(c,R-c,y,D-BQ)P^ +(AQ+o,T +CD-BR-o,y,U)P” 
+(AR+CU - Moy, -BT)P' +CM + AT =0 
On simplification we have the four degree polynomial function in terms of P` given by 
a, P” +a P” +a,P” +a,P* +a, =0 (29) 
where 
a, =0,Q 
a, =(o,R-o,7,D-BQ) 
a, =(AQ+o,T+CD-BR-o,y,U) 
a, =(AR+CU - Mo,y, — BT) 
a, =CM + AT 
Substituting P' from equation (17) into (29) we have the expression for the endemic 


equilibria which satisfy the following four degree polynomial in terms of Р, given by 
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YF; +Y F +Ү,Е; +Ү,Е; +Y, =0 (30) 
where 
Y, = ayrz + ar? + jr; + Gyr, +a, 
Y= lar? +2a,rê +3a,r, +4a, k 
Y, = (a.r? +3a,r„k? +6a, k? 
Y, = (азгь +4a,) 2 
Y, =a,k^ 
Since the variables $5,,1,,,1,,,S,. li Lp F, are expressed in terms of P then we 
substitute P` =r, Fi/(k+F;) in equation (29) we obtain the polynomial function with 
degree four in (30) which represents the presence of endemic equilibrium points. This 
shows that there are possible four roots for F; which further implies that there are at most 


four possible endemic equilibrium points. 


4.5 Global stability of endemic equilibrium point 

Global stability of the endemic equilibrium (EE) is investigated using Lyapunov method 
and LaSalle's invariance principle. The approach which has been found to be useful for 
compartmental epidemic models with any number of compartments (Korobeinikov and 
Maini, 2004). To achieve this we construct a suitable Lyapunov function of the form: 


8 
L-Y A(x; - x Inx,) (31) 
i=l 
where A, is a properly selected constant, x, is the population of i” compartment, x, is the 


equilibrium value of x, and А » 0. The Lyapunov function denoted by L is continuous 


and differentiable. We have: 


Нара Е =A (85 58 Ing; + A (Ij Fells} 
+ A (Ly, Тш, )+ AS, = 5а, 
PAN = Dll БАГ) 
+ A (Le -LplnL, )+ A (Fs — Еа, ) 


(32) 


The global stability of endemic equilibrium (EE) holds if its time derivative T <0. 


Proof: 
The time derivative of the Lyapunov function L is given by 


By adopting the approach by McCluskey, (2006) we have the following expressions 


-A- 
+A 
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7(Ф) is non-positive by following the approach of Korobeinikov (2002, 2004, 2007) and 


McCluskey (2006). Thus 7(Ф)<0 for all Ф>0 Hence, 4-<0 іп Ф and is zero 
t 


when Ф = Ф“. Therefore the largest invariant set in Ф such that 22 = 0 is the singleton 


(Ф )which is our endemic equilibrium point. By LaSalle’s invariant principle (LaSalle, 


1976) we conclude that (D`) is globally asymptotically stable (GAS). Thus, we establish 
the following Theorem 
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Theorem 5: The endemic equilibrium point D” of model system (3) is globally 
asymptotically stable in Ф if R, >1 and unstable otherwise. 


5 Simulations on the basic model state variables over time 

The simulation of the model is conducted to find out the dynamics of Tungiasis disease in 
the population when there is no intervention. It is performed using MALAB, and we set 
time in weeks and days. We use parameter values whose sources are from literature and 
others are estimated. Furthermore through numerical simulation we illustrate the stability 
of the endemic equilibrium states in human, animal reservoirs and flea in the environment 
whenever they exist. The parameter values used for simulation are shown in Table 3. 


5.1 Parameter Estimation 
In this section we estimate some of the parameter values of Tungiasis model. We estimate 
the human mortality rate w,,, taking into consideration that the average life expectancy of 


the human population in Tanzania is 60.9 years (UNICEF, 2015) which is equivalently to 
4.5x10? per day. According to Tanzania population (TP, 2016) the human birth rate is 
0.00011 per day. The maturation rate y, of larvae into adult flea is estimated to be 0.0105 


рег day and the shedding rate of fleas е, to be 0.4 per day. The death of the flea occurs 


around day 25 post-penetration (Eisele et a/, 2003), we therefore assume the life span of a 
flea to be 25 days, which implies that its death rate >, is 0.04 per day. The concentration of 


jigger fleas in the soil environment is not known; we therefore consider the number of fleas 
in one cubic meter of sand and assume the maximal larvae carrying capacity K to be 
1x10?cell/m? and the half saturation constant for fleas k to be1x10* сеП/ т’. The disease 
induced death rate for animal reservoir population с; is estimated to be 0.037 per day and 
the disease induced death rate for human population o,, is estimated to be 0.011 per day. 
The natural death rate of animal reservoir ш, ranges from (360—3600)days (Radostits, 
2001; Gaff et al, 2007). The values of effective contact rates are p, =0.26 per day 
(Allerson et al, 2013), р, =0.052 perday (Gaff et al, 2007), £,,,=0.19perday and 
Pra = 0.48 per day are estimated. 
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Table 3: Parameter values used for the simulation of model variables 










































































Parameter Value/Range Source/References 
k 1х10* cell/m? Estimated 
K 1x10? cell/m* Estimated 
V ï 0.0105per day Estimated 
Он 0.011 per day Estimated 
с, 0.037 per day Estimated 
Hy 0.000045 per day UNICEF. (2015) 
H 0.0028(360— 3600) ' per Gaff et al. (2007); Radostits. 
day (2001) 
p» 0.04 per day Eisele et al (2003), 
Fp 0.58 per day Estimated 
L 0.08 per day Estimated 
Шар 0.19 рег day Estimated 
Bea 0.48 per day Estimated 
Pau 0.052 per day Gaff et al. (2007) 
Pa 0.26 (0.091-0.9) per day Allerson et a/. (2013) 
by 0.00011 per day TP, (2016). 
b, 0.022 per day Gaff et al. (2007) 
£4 0.40 per day Estimated 
ó, 0.12 per day Estimated 
Gi 0.4 Estimated 
Ara 0.6 Estimated 
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5.2 Discussion of the results 
Figure 2 illustrate the dynamics in humans, animal reservoirs and flea populations, 
showing the behavior of eight state variables S$,,7,,,1,,, $,,1,,1,,, Г апа F; over time 


when there is no intervention. 
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Figure 2: Dynamics in model sub-populations 


From Figure 2: The red dotted line is the larvae population and the blue dotted line is the 
adult flea population. The larvae population decreases exponentially to attain endemic 
equilibrium point as they die naturally and as they mature into adult fleas. The population 
of adult fleas generated from larvae maturation increase exponentially to a certain point in 
time then decrease to attain endemic equilibrium point as they die a natural death and as 
they burrow into the epidermis of the hosts. The red line is the susceptible humans, the 
black line is mildly infested humans and the cyan line is the severely infested humans. The 
susceptible humans decrease exponentially due to natural death and as they acquire 
infestation from severely infested animal or from the flea infested environment and 
eventually attain the endemic equilibrium point. The decrease in susceptible human 
population results into the growth of the mildly infested humans that grow exponentially 
after sometime they start declining due to natural death and as they acquire more 
infestation from the environment to join the severely infested human class then eventually 
attains endemic equilibrium point. The severely infested humans grow exponentially to a 
certain point in time then start declining due to natural and disease induced death and 
finally attains endemic equilibrium point. The blue line is the susceptible animal reservoirs, 
the green line is the mildly infested animal reservoirs and the magenta line is the severely 
infested animal reservoirs. The susceptible animal reservoirs decrease exponentially due to 
natural death and as they acquire infestation from severely infested animal reservoirs or 
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from the flea infested environment. The decrease of susceptible animal reservoirs results 
into the growth of the mildly infested animal reservoirs that grow exponentially and then 
start declining due to natural death and as the acquire more infestation from the 
environment to join severely infested animal reservoirs class then eventually attain 
endemic equilibrium point. The severely infested animal reservoirs grow exponentially to a 
certain point in time then start declining due to natural and disease induced mortalities 
and finally they attain endemic equilibrium point. 


5.3 The variation of population variables on the dynamics of Tungiasis over time 
The numerical simulations as illustrated in Figures 3(a), 3(b), 3(c), 4(a), 4(b), 4(c), 5(a) 
and 5(b) show the variations of the model variables 8,,1,,1,,,9,,1,,1,,, Lg and F; over 


time. The trajectories of the model variables originating from different initial values 
converges to a common endemic equilibrium levels which implies the existence and 
stability of the endemic equilibrium states in human, animal reservoirs and flea 
populations whenever they exist. 
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Figure 3(a): Stability of the endemic equilibrium for susceptible human population 
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Figure 3(b): Stability of endemic equilibrium for mildly infested human population 
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Figure 3(c): Stability of the endemic equilibrium for severely infested human population 
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Figure 4(a): Stability of the endemic equilibrium for susceptible animal population 
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Figure 4(b): Stability of the endemic equilibrium for mildly infested animal population 
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Figure 4(c): Stability of the endemic equilibrium for severely infested animal population 
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Figure 5(a): Stability of the endemic equilibrium for jigger larvae population. 
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Figure 5(b): Stability of the endemic equilibrium for adult jigger flea population 
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6 Conclusions 

The Tungiasis model with sub-models of humans, animal reservoirs and flea populations, 
has been developed and analyzed to study the tungiasis dynamics in endemic areas. The 
basic reproduction number R, computed was used to determine the stability of the disease- 


free equilibrium point (DFE). We then proved the existence and stability of endemic 
equilibrium point (EE) using Lyapunov direct Method combined with LaSalle’s Invariance 
Principle. The analytical results show that the DFE is locally asymptotically stable, if R, <1 


and the EE is globally asymptotically stable when R, >1. From numerical simulations we 


have observed that without intervention the populations vary for some time but ultimately 
approach the endemic equilibrium levels in the long run which imply the existence and 
stability of endemic equilibrium point. 

This work provides a basic dynamical model that can be used to understand the 
transmission dynamics of tungiasis transmission that involves the interaction between 
human population, animal reservoir population and the flea infested environment. 
Understanding tungiasis transmission dynamics will help to design proper control 
strategies and evaluate their potential impact in reducing tungiasis morbidity and 
mortality in the endemic communities at large. 
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